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Abstract. We present the Josephson junction intersected superconducting trans- 
mission line resonator. In contrast to the Josephson parametric amplifier, Josephson 
bifurcation amplifier and Josephson parametric converter we consider the regime of 
few microwave photons. We review the derivation of cigenmodc frequencies and zero 
point fluctuations of the nonlinear transmission line resonator and the derivation of 
the eigenmode Kerr nonlinearities. Remarkably these nonlinearities can reach values 
comparable to Transmon qubits rendering the device ideal for accessing the strongly 
correlated regime. This is particularly interesting for investigation of quantum many- 
body dynamics of interacting particles under the influence of drive and dissipation. 
We provide current profiles for the device modes and investigate the coupling between 
resonators in a network of nonlinear transmission line resonators. 
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1. Introduction 

The theory of interacting quantum many-body systems is a vibrant discipline in 
physics full of unsolved riddles like high temperature superconductivity and promising 
technological prospects such as topologically protected quantum states, novel schemes 
for quantum error correction or one way quantum computing. Minimal Hamiltonians 
are considered for both situations: to explain quantitatively phenomena in nature or to 
provide prescriptions how to build artificial systems for technological applications. 

There exists a constantly growing amount of analytical and numerical tools to 
investigate the properties of many body Hamiltonians. However the analytical tools 
typically provide answers only in limited regimes while the numeric approaches are 
constrained by the exponential growth of the Hilbert space with the number of particles. 
Dating back to an idea of Feynman pQ who proposed to use an computer that uses the 
principles of quantum mechanics itself to investigate quantum mechanical many-body 
Hamiltonians the idea of quantum simulation [2J is gaining a lot of attention nowadays. 
A quantum simulator is a well controllable quantum systems that emulates the physics 
of other systems which are less amenable to experimental investigation. Therefore 
quantum simulators are typically strongly enlarged versions of the "real world" system 
we try to investigate. This provides us with the possibility of tuning parameters of 
the Hamiltonian and with individual addressability for measurements because of larger 
lattice constants. Quantum many-body Hamiltonians can for example be simulated 
with cold atoms trapped by laser fields in various shapes and dimensions [3j, in ion 
traps [U |5] or arrays of cavity quantum electrodynamics (QED) systems [HI [7]. 

Massless photons in cavities don't interact. However by strongly coupling the 
photons to nonlinear systems, joint excitations emerge, called polaritons [HI [7J El l9l [TO], 
and the nonlinearity of the onsite spectrum for the polaritons can be interpreted 
as an interaction (TTJ H2] . Cavities can be coupled in arbitrary topologies thereby 
generating networks for moving photons. The polaritons inherit this abiltiy from the 
photons. The description of the fundamental excitations of the quantum simulator 
in terms of polaritons is only valid in the strong coupling regime where the coupling 
between the cavity and the nonlinear system exceeds the decay rates of both. In circuit 
quantum electrodynamics (circuit QED), the microwave realization of cavity quantum 
electrodynamics, one can easily achieve the strong coupling regime. In circuit QED 
superconducting wires (transmission line resonators) that confine microwave photons 
are coupled to Josephson junction based lumped element circuits ( Josephson qubits) [13, 
HU [T3]. Strong zero point fluctuations of the field due to the restricted dimensionality 
of the transmission line resonator and the artificial coupling to the Josephson qubit 
result in exceptionally high coupling rates [TH] while the superconducting gap ensures 
low dissipation. The polariton approach has two major flaws. Firstly, the ratio of 
nonlinearity to hopping strength is changed by changing the polariton from purely 
photonic to Josephson qubit excitation. Therefore nonlinearity and hopping are 
mutually exclusive. Secondly, we always have two polariton species separated in energy 
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space by the coupling strength which provides us with unwanted but unavoidable 
crosstalk. 

In this work, we propose making the resonator itself nonlinear by intersecting it with 
a Josephson junction [T7J [HI [19]. We show that the modes of such a device can achieve 
amounts of nonlinearity comparable to transmon qubits [201 [21] for realistic circuit 
parameters and that the mode separation in frequency space is sufficient to suppress 
the crosstalk that is introduced by the nonlinearity. Moreover the coupling of these 
resonators is independent of the nonlinearity. Notably, a set of experimental analysis 
methods for the quantum statistics of resonator networks have been demonstrated in 
recent years [221 1221 121] • These properties make the nonlinear transmission line resonator 
an ideal building block for a Bose Hubbard quantum simulator. 

The article is organized as follows. First we review the derivation of frequencies and 
zero point fluctuations of the eigenmodes of the nonlinear transmission line resonator 
[18] . Then we derive the individual Kerr nonlinearities for the eigenmodes in a low energy 
and rotating wave approximation. Finally we provide the current profiles including 
current through the inductive part of the Josephson junction and current through the 
transmission line resonator and discuss capacitive coupling of the eigenmodes between 
adjacent transmission line resonators. 



2. Nonlinear transmission line resonator 



We consider a transmission line resonator intersected in the middle by a Josephson 
junction. In this section we briefly review the decomposition of the microscopic 
Lagrangian of the nonlinear transmission line resonator into independent modes, we 
provide some physical intuition for the modes and calculate their coupling strength in 
a network of nonlinear transmission line resonators. 



2.1. spectrum 

The microscopic Lagrangian of the whole setup reads, 

£ = £ l + £ r + £j J , (1) 

with the Lagrangian of the transmission line to the left and the right of the Josephson 
junction, 

o/f 

A/r= / l^l/rf ~ ^{d^/rfdx , (2) 
-I/O 

with the flux field (f){x) = V(x,t')dt' [22] and the capacitance and inductance per 
unit length of transmission line resonator c and /. The Lagrangian of the Josephson 
junction £jj depends on the flux drop at the Josephson junction S(f> = (f>i\ x =o — (fi r \x=o, 

£jj=^fS<j> 2 + Ej cos A<f>), (3) 
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where Cj and Ej are the capacitance and Josephson energy of the junction respectively. 
0o = h/2e is the quantum of flux. We proceed by separating the linear from the purely 
nonlinear parts of the Lagrangian. This way we can exactly diagonalize the linear parts 
and get the normal modes of the system. We then reintroduce the nonlinear parts as an 
perturbation to the oscillator modes of the system. It suffices to consider the nonlinear 
cosine inductance of the Josephson junction when separating the nonlinear parts from 
the linear ones , 



^<^-^ 2 H 1 + gS ( >1' (4) 



where Lj = (^/{At^Ej) is the Josephson inductance. Considering only the linear part 
of the Lagrangian, the Euler Lagrange equations provide us with the wave equation in 
the two arms of the transmission line resonator, 

9t<Pl/r = V 2 d 2 x <pl/ r , (5) 

where the phase velocity is given by the inductance and capacitance per unit length 
v = l/yfd. Boundary conditions for the flux field can either be derived by Kirchhoffs 
rules or also via the Euler Lagrange equations. The current in the transmission line 
resonator is proportional to the spatial derivative of the flux field / = jd x (p and has to 
vanish at the ends of the transmission line, 

d x <pi/ r \ x=T L = 0. (6) 

We also know that the current flowing into the Josephson junction from one side has to 
exit on the opposite side, 

d x (pi\ x =o = d x 4> r \ x= o . (7) 

The current flowing through the junction can either flow into the shunting capacitor or 
the Josephson junction. The current flowing through the junction is related to the flux 
drop across the junction in the Josephson constitutive relation lj = I c sm((2n /(f> )S(f)). 
The current into the capacitor in turn reads, I cap = Cj5(j). For the linearized Lagrangian 
the condition for the current flowing through the Josephson junction reads, 

-\d x <f> l = Cj8'^+^8(l>. (8) 

To find the spatial eigenmodes of the nonlinear transmission line resonator we first make 
an ansatz with separation of variables, <j)(x, t) = f(x)g(t) with fi(x) = A t cos(k(x+L/2)) 
for the transmission line to the left of the Josephson junction and f r (x) = A r cos(k(x — 
L/2)) for the right part of the transmission line. With this ansatz the necessity of 
vanishing currents at the ends of the transmission line is automatically satisfied and the 
wave equation is fulfilled provided that oo = kv , where u is the frequency of the temporal 
part of the flux function g + u 2 g = 0. There are two distinct possibilities to fulfill the 
boundary conditions at the Josephson junction: we could either choose Ai = A r which 
provides us with the symmetric modes of the resonator or we could choose A\ = —A r to 
get the antisymmetric modes. The symmetric modes show no flux drop at the Josephson 
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Table 1. 

transmission line 

characteristic impedance Z$ 50 £1 

phase velocity v 0.94-10 8 m/s 

inductance per length I 5 ■ 10~ 7 H/m 

capacitance per length c 2-10~ 10 F/m 

Josephson Junction 

shunting capacitance Cj 1.9 ■ 10~ 12 F 



junction and consequently there is no current flowing through the junction. Therefore 
the symmetric modes of the transmission line resonator are completely unaffected by 
the presence of the Josephson junction and we omit them in the further discussion of 
the system dynamics. For the antisymmetric modes we choose the amplitudes to be 
A\ = —A r = 1 and get the following transcendental equation for the mode frequencies 



,2 



where u p = 1/y/LjCj is the Josephson plasma frequency. The solution of the above 
transcendental equation is the spectrum of the system, not taking into account the 
symmetric modes, presented in figure [I] a) for generic parameters given in table [l] We 
plotted the spectrum as a function of the critical current I c of the Josephson junction. 
The resonator modes without the Josephson junction do not depend on the cricital 
current of the Josephson junction but the plasma frequency cu p = yTji^oCj) does. 
If the plasma frequency and the frequency of the antisymmetric modes of the bare 
transmission line resonator (without Josephson junction) are detuned the frequencies 
of the combined system are the frequencies of its constituents, Josephson junction or 
transmission line resonator. But if the plasma frequency is resonant with one of the 
symmetric modes we get an anticrossing of the two frequencies which is of the order of 
the mode frequency itself accounting for the large coupling between Josephson junction 
and transmission line resonator. 

We deduce the following generalized "scalar product" between the spatial modes 
with the help of the transcendental equation of the mode frequencies, 

c / f n {x)f m (x)dx + C s 5f n 5f m = 5 niin r] n , (10) 



with, 



>!,,='■ I i + S S^(i + $;) ) (ID 



2 kl 2Lj V ^1 



where 5f n = f n ,i\x=o ~ fn,r\x=o is the normalized flux drop of the spatial modes. We 
use this "scalar product" to decompose the linearized Lagrangian into independent 
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Figure 1. a) Numerically calculated frequencies of the linear part of the transmission 
line resonator intersected by a Josephson junction. Anti-crossings show up where the 
plasma frequency of the Josephson junction lo p (dashed ascending graph) matches 
one of the frequencies of the odd modes of the transmission line resonator (dashed 
horizontal graphs), b) Nonlinearity parameter of the different modes U n . If the 
frequency of the mode is near the plasma frequency of the Josephson junction, 
the nonlinearity e 2 /(4C s ) is inherited, c) Coupling between different antisymmetric 
modes of adjacent nonlinear transmission line resonators. Comparison with the 
mode frequencies reveals the connection of coupling strength g n = (C c w„)/(4?7„) and 
frequency of the respective mode. Dashed vertical lines mark the values of the critical 
current I c where we computed the current profile of the mode cf. figure [2] 
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oscillators of frequencies u n and masses i] n . After a Legendre transformation and 
omitting constant terms we get the Hamiltonian, 



2n 



Up to here everything we did was classical physics, now we proceed by quantizing the 
Hamiltonian in the usual way by introducing independent raising and lowering operators 
for the different modes, 



; ~ l\ I — Y~ ( a « ~~ a n. 

1 / u 

9n = (On + <) , 

y/2r] n UJ n 



with [a n , aJJ = 5 n ,m- Finally, we arrive at the following Hamilton operator, 



n=l 



with 



H' M)ulm = <o>(V>) • yvr) . (12) 



S(f> = ^2\ n (a n + a\) and, A n 



The nonlinearity again reintroduces interactions between the modes and nonlinearities 
for the individual modes. As we will show below these are all much smaller than the 
differences between the mode frequencies. They only play a role in strong driving 
scenarios like the Josephson parametric converter. We are only interested in few 
microwave photon physics and therefore may omit the residual coupling in a rotating 
wave approximation. The individual mode nonlinearities however can not be omitted. 
For their derivation we start to decompose the cosine term with a product-to-sum 
formula which is applicable in this case because all individual flux operators commute, 

(oo \ / oo 

X n (a n + 4) J = cos (Ai(gh + a})) cos j \ n (a n + a n ) 
n=l / \n=2 

(oo 
^ A n (a n + al) 
n=2 

The sine terms only contain rotating operator products and therefore can be neglected 
in the framework of the rotating wave approximation. If we further proceed this way 
we may transform the cosine of the sum of all mode flux operators into the product of 
the cosines of all mode flux operators, 



cos I y X n(a n + 4) J ->■ Yl cos (A n (a n + ajj) 

\ n / n 
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The cosines of the mode flux operators may be rewritten in a sum of normal ordered 
products of mode operators. We may further also omit all operator products containing 
a different amount of raising and lowering operators, again because of the rotating wave 
approximation. Finally we get as a rotating wave, low energy approximation for the 
cosine of a mode flux operator, 



cos 



A 4 

(A n (a n + a]J) -> ( 1 - \ 2 n a ] n a n + -^a\a\a n a n + ... 



We concentrate on the fundamental mode and consider all the other modes to be in 
the vacuum state which leads us to an Hamiltonian only describing the fundamental 
antisymmetric mode of the transmission line resonator coupled to the plasma mode of 
the Josephson junction, 

Hi = (cji — 5u) a\a\ — U a\a\a\ai , 



with 



5u = EjXl 1 - J] 




e 

78=1 

/ / [[< . 

71=1 

where 5u is a small renormalization of the fundamental mode frequency due to the 
nonlinearity. The Kerr nonlinearity U is plotted in figure [l] b) as a function of the 
Josephson junctions critical current I c . Direct comparison with the spectrum of the 
system reveals the origin of the nonlinearity. For values of I c where the respective 
frequency is defined by the plasma frequency of the Josephson junction, the mode 
inherits the full nonlinearity of the Josephson junction. 



2.2. modes 

We now decomposed the linear part of the Lagrangian into independent harmonic 
oscillators with frequencies u n and mode capacitances rj n . The latter only make sense in 
combination with the normalization convention we chose for the spatial mode functions 
f n (x). In order to get the correct physical intuition of the oscillation mode, we have to 
revert to the well known quantities of oscillating circuits, currents and charges or rather 
coulomb potential. We start by calculating the current through the Josephson junction. 
The observable for the current through the Josephson junction is, 

ij = I e em(8$). (13) 

Lets suppose that all modes except for the fundamental mode are in their vacuum state 
and the fundamental mode is in a Fock state \n0 . . . 0). In this state the mean value of 
current flowing through the Josephson junction vanishes comparable to the vanishing 
displacement of a quantum harmonic oscillator in a Fock state. The variance of the 
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Figure 2. Variance of the current in the transmission line resonator and Josephson 
junction for the first Fock state of the fundamental mode. 



current however doesn't vanish and provides us with an estimation of the mean current 
flowing through the Josephson junction, 



Aij = 7 c y (nO . . . 0| sm(Scj)) 2 \n0 . . . 0} . (14) 

With the same arguments as presented above for the derivation of the Kerr parameter, 
we calculate the variance, which is exact up to the second Fock state, 



AI, 



\ - n e ~ 2xi ^ - 4A ? n + 4A ' n ( n - • 



(15) 



The observable for the current in the transmission line to the left and the right of the 
Josephson junction is, 



i r (x) = jd x (j)(x) = 



d x f n {x) 



(16) 



Here again the mean current in the transmission line resonator for a Fock state vanishes 
and we evaluate the variance of the current, 



AI r 



I 



\d x fi{x)\ 2 + yl f \d x f k (x)\ 2 



(17) 



figure [2] shows the variance of the current in the left half of the nonlinear transmission 
line resonator and the inductive part of the Josephson junction. The variance of the 
current in the right part is for symmetry reasons the mirror image of its counterpart in 
the left half. 
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2. 3. coupling 

Next we are calculate the strength of the coupling between two neighboring transmission 
line resonators in a network. We only consider capacitive coupling of resonators where 
the respective ends of the central lines of adjacent transmission line resonators are either 
close to each other or connected by interdigitated capacitors. To integrate the coupling 
into our theoretical model we have to include the energy of the coupling capacitor, 
with coupling capacitance C c , into the Lagrangian of the nonlinear transmission line 
resonator, 




where 0i/2 is the flux field of adjacent nonlinear transmission line resonators. After a 
Legendre transformation to get the corresponding energy term in the Hamiltonian we 
would get different conjugate momenta Tr n . We neglect this effect because the coupling 
capacitance C c is very small compared to the overall capacitance of the nonlinear 
transmission line resonator and get 0i/ 2 = —i J2 n ■\/u n /{2'r) n )f n {x) (ai/ 2 , n — a \/2,n) • 
The coupling term in the Lagrangian therefore provides us with many different effects: 
We get renormalizations of the different resonator mode frequencies, exchange couplings 
between different modes of adjacent resonators and on the same resonator, which we 
neglect in a rotating wave approximation, and exchange coupling of the same mode of 
adjacent resonators, 

Hn,c = -— ^— (a\ n a 2 ,n + «l,n4n) ■ ( 19 ) 
4 T] n \ ' 'J 

The dimensionless capacitance of the fundamental modes of the nonlinear transmission 
line resonator is not affected by a change of the Josephson junction's critical current 
for generic parameters (cf. table [T]) of our setup . The coupling g„ = (C c co n ) /(4r/ n ) of 
the same modes in adjacent resonators is therefore determined by the frequency of the 
modes cf. figure [I] c), where we plotted the coupling of the first three antisymmetric 
modes of the nonlinear transmission line resonator. This in turn enables us to increase 
the nonlinearity while increasing the coupling, provided we do not use the fundamental 
mode of the nonlinear transmission line resonator. 



3. Summary 

We reviewed the derivation of the eigenmode frequencies and respective zero point 
fluctuations for the Josephson junction intersected transmission line resonator. We 
derived the Kerr nonlinearities of the individual eigenmodes and calculated the current 
in the transmission line resonator and the Josephson junction. Finally we showed 
that the coupling depends on the eigemode frequency and is thereby independent of 
the nonlinearity. This lack of mutual exclusivity of nonlinearity and coupling makes 
quantum simulators with nonlinear transmission line resonators superior to polariton 
approaches. 
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